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ABSTRACT 

We have made the first measurement of the double- inclusive B/B energy distribution 
in e + e~ annihilations, using a sample of 400,000 hadronic Z°-decay events recorded in 
the SLD experiment at SLAC between 1996 and 1998. The small and stable SLC beam 
spot and the CCD-based vertex detector were used to reconstruct -B/5-decay vertices 
with high efficiency and purity, and to provide precise measurements of the kinematic 
quantities used to calculate the B energies in this novel technique. We measured the 
B/B energies with good efficiency and resolution over the full kinematic range. We 
measured moments of the scaled energies of the B and B hadrons vs. the opening angle 
between them. By comparing these results with perturbative QCD predictions we tested 
the ansatz of factorisation in heavy-quark production. A recent next-to-leading order 
calculation reproduces the data. 



* Work supported in part by Department of Energy contract DE-AC03-76SF00515. 



1 Introduction 



The production of heavy hadrons (H) in e + e~ annihilation provides a laboratory for the 
study of heavy-quark (Q) jet fragmentation. This is commonly characterized in terms 
of the observable xjj = 2Eh / \fs, where Eh is the energy of a B or D hadron containing 
a b or c quark, respectively, and y/s is the cm. energy. The distribution of xh, D(xh), 
is conventionally referred to as the heavy-quark 'fragmentation functionfl 

In recent publications we presented [l], [| the results of a new method for recon- 
structing .B-hadron decays, and the B energy, inclusively, using only charged tracks, 
in the SLD experiment at SLAC. We used the upgraded charge-coupled device (CCD) 
vertex detector, installed in 1996, to reconstruct S-decay vertices in Z° decays with high 
efficiency and purity. Combined with the micron-sized SLC interaction point (IP), our 
precise vertexing allowed us to reconstruct the total transverse momentum of the tracks 
from 5-decays, and therefore the transverse momentum and mass associated with the 
neutral particles in the -B-decays. This allowed us to reconstruct accurately the energy 
of B hadrons. These studies yielded the most precise measurement of the 6-quark frag- 
mentation function, and allowed us to test models of heavy-quark fragmentation. Of 
the 9 models tested, only 4 were consistent with our precision data at better than the 
1% level based on a x 2 probability. This allowed us to reduce the model-dependent 
systematic uncertainty on the 6-quark fragmentation function. 

We have extended these studies and applied similar 'topological' vertexing techniques 
to tag events in which we reconstructed the energies of both leading B hadrons produced 
via e + e _ — > bb — > BB + X. We measured the moments of the single-inclusive 5-hadron 
scaled-energy distribution dN/dxs'- 

as well as the moments of the double-inclusive scaled-energy distribution: 

Aj(0) = / / x % bIx^-^- 7 dx m dx B 2, (2) 

J J N d dx Bidx B2d cos <p 

where xbi and xb2 are the scaled energies of the two B hadrons and the label is arbitrary, 
<p is the angle between their flight directions, and i and j are integers >1. We formed 
the normalised moments: 

G l3 {4>) = IhjiOl il);!);) (3) 

Following the method proposed in 0] we used these quantities to test the ansatz of 
factorisation as applied to perturbative Quantum Chromodynamics (pQCD) calculations 
of e + e~ — > bb events. 

* Unless stated otherwise, in these studies we do not distinguish between hadrons and antihadrons. 
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Due to 'soft', or non-perturbative, effects, Di and Dij cannot be predicted absolutely 
in pQCD. Rather, the respective pQCD calculation must be folded with models of the 
non-perturbative ('hadronisation') process in order to derive predictions that can be 
compared with experimental data. However, provided that the ansatz of factorisation 
holds ||, namely that calculation of the perturbative and non-perturbative phases can 
be separated by (an arbitrary) factorisation scale //p, the dependence on //p cancels in 
Eq. H and hence G^ can be calculated absolutely in pQCD, up to possible 'higher twist' 
effects of order 1 / ■sfs, with no dependence on hadronisation models. Comparison of the 
measured Gij with pQCD predictions hence allows both a test of this ansatz and of the 
perturbative calculations. The can be derived at leading order (LO) in pQCD using 
the numerical results in ||. In addition, next-to-leading order (NLO) predictions for G^ 
have been calculated recently ||]. 

In Section 2 we describe the detector and the selection of e + e~ — ► hadrons events used 
in this analysis. We present in Section 3 the first measurement of the double-inclusive 
5-hadron energy distribution and, in Section 4, of the normalised moments Eq. ^. In 
Section 5 we describe the estimation of the errors on our measurements. Finally, in 
Section 6, we test the ansatz of factorisation and compare our data with the pQCD 
predictions. 

2 Apparatus and Hadronic Event Selection 

This analysis is based on roughly 400,000 hadronic Z° events produced in e + e~ an- 
nihilations at a mean center-of-mass energy of y/s = 91.28 GeV at the SLAC Linear 
Collider (SLC), and recorded in the SLC Large Detector (SLD) between 1996 and 1998. 
A general description of the SLD can be found elsewhere ||. The trigger and initial 
selection criteria for hadronic Z° decays are described in Ref. ||. This analysis used 
charged tracks measured in the Central Drift Chamber (CDC) ]7j and in the upgraded 
Vertex Detector (VXD3) ||. Momentum measurement was enabled by a uniform axial 
magnetic field of 0.6T. The CDC and VXD3 gave a momentum resolution of <J p± /p± = 
0.01 ©0.0026p_i_, where p± is the track momentum transverse to the beam axis in GeV/c. 
In the plane normal to the beamline the centroid of the micron-sized SLC IP was recon- 
structed from tracks in sets of approximately thirty sequential hadronic Z° decays to a 
precision of a r jp ~ 4 /im. The IP position along the beam axis was determined event by 
event, using charged tracks, with a resolution of a z IP ~ 20 /im. Including the uncertainty 
on the IP position, the resolution on the charged-track impact parameter (d) projected 
in the plane perpendicular to the beamline was a 1 ^ = 8©33/(psin 3//2 9) //m, and the 
resolution in the plane containing the beam axis was o z d = 10©33/(psin 3//2 9) //m, where 
9 is the track polar angle with respect to the beamline. 

A set of cuts was applied to the data to select well-measured tracks and events 
well contained within the detector acceptance. Charged tracks were required to have a 
distance of closest approach transverse to the beam axis within 5 cm, and within 10 cm 
along the axis from the measured IP, as well as | cos6>| < 0.87, and p± > 0.15 GeV/c. 
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Events were required to have a minimum of five such tracks, a thrust axis || polar angle 
w.r.t. the beamline, 6?, within | cos6t\ < 0.80, and a charged visible energy E vis of at 
least 20 GeV, which was calculated from the selected tracks assigned the charged pion 
mass. The efficiency for selecting a well- contained Z° — > qq(g) event was estimated to be 
above 97% independent of quark flavor. The selected sample comprised 313,447 events, 
with an estimated 0.10 ± 0.05% background contribution dominated by Z° — > t + t~ 
events. 

For the purpose of estimating the efficiency and purity of the selection procedures 
we made use of a detailed Monte Carlo (MC) simulation of the detector. The JETSET 
7.4 |L0] event generator was used, with parameter values tuned to hadronic e + e~ annihilation 



data [11 1, combined with a simulation of I?-hadron decays tuned [12] to T(45) data 
and a simulation of the SLD based on GEANT 3.21 |I3"| . Inclusive distributions of 
single-particle and event-topology observables in hadronic events were found to be well 
described by the simulation 0. Uncertainties in the simulation were taken into account 
in the systematic errors (Section ||). 



3 5-Hadron Selection and Energy Measurement 

The event sample for this analysis was selected using a 'topological' vertexing technique 
based on the detection and measurement of charged tracks, which is described in detail in 
Refs. |14|, [15], [16|] . We considered events in which we found decay vertices corresponding 



to both the leading B and B hadrons. First, the Durham algorithm [17] was applied to 
the selected hadronic events, with a y c parameter value of 0.015, in order to define a jet 
structure in each event. We found that this algorithm and y c value minimized the number 
of B (and D) decay tracks assigned to the wrong jet. This is an important feature for 
our analysis because we used vertex-related variables derived only from charged tracks. 
Events containing 2, 3, or 4 jets were retained for further analysis. 

In each selected event, the vertexing algorithm was applied to the set of tracks in 
each jet. Vertices consistent with photon conversions or K° or A decays were discarded. 
Events were retained in which a vertex was found in exactly two jets. 35,137 events were 
selected, of which 89.4% were estimated to be of bb origin. The efficiency for selecting 
true bb events was estimated to 36.3%. 

The large masses of the B hadrons relative to light-flavor hadrons make it possible to 
distinguish S-hadron decay vertices from those vertices found in events of lighter flavors 
using the vertex invariant mass, M. However, due to those particles missed from the 
vertex, which are mainly neutrals, M cannot be fully determined. M can be written 

M = sjKh + P? + PcM + V^o + p t + P~oi (4) 

where M c h and M are the total invariant masses of the set of vertex-associated tracks 
and the set of missing particles, respectively. P t is the total charged-track momentum 
transverse to the B flight direction, which, by momentum conservation, is identical to 
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the transverse momentum of the set of missing particles. P c m and Pqi are the respective 
momenta along the B flight direction, which we take to be the vector joining the IP to 
the vertex. 

The lower bound for the mass of the decaying hadron, the 'Pt-corrected vertex 
mass' PI, 



M Pt = y/M& + P?+\P t \ (5) 

was used as the variable for selecting B hadrons. Figure [l] shows the distribution of Mp t 
for vertices in the selected event sample, and the corresponding simulated distribution. 
Events were selected that contained at least one vertex with Mp t > 2.0 GeV/c 2 and 
Mp t < 2 x M c h. The latter cut was found to reduce the contamination from fake 
vertices in light-quark events. 

In order to improve the bb purity of the sample, events were selected in which both 
vertices had a flight length, d vtx , such that 0.1 < d vtx < 2.3 cm; in which at least one 
vertex contained two 'significant' tracks, i.e. tracks with a normalised impact-parameter 
significance, d/<Jd, of at least 2 units; and in which the angle between the vertex flight 
vectors, 0, satisfied cos0 < 0.99. The last cut was effective at removing events in which 
either a gluon had split into heavy quarks, or the jet-finder had artificially split a single 
heavy-quark jet into two jets. 

The energy of each B hadron, Eb, can be expressed as the sum of the reconstructed- 
vertex energy, E c h, and the energy of those particles not associated with the vertex, E . 
We can write 

El = Ml + P 2 + P 2 (6) 

The two unknowns, M and Pqi, must be found in order to obtain E Q . One kinematic 
constraint can be obtained by imposing the 5-hadron mass on the vertex, Mj, = E 2 B -Pl, 
where Pp = P c m + Pqi is the total momentum of the B hadron. From Equation (B) we 
derive the following inequality, 



M 2 ch + P t 2 + VM 2 + P? < M B , (7) 

where equality holds in the limit where both Pqi and P c u vanish in the B hadron rest 
frame. Equation (|7|) effectively sets an upper bound on M , and a lower bound is given 
by zero: 

< Ml < Ml max , (8) 

where 



M 2 0max = Ml - 2Mb yM 2 h + P 2 + M 2 h . (9) 

Because M peaks near M 0max , fl5j we set Ml = Ml max if Ml max >0, and Ml = 
if Ml max <0. We calculated P ol : 



Ml - ( Mj + P 2 ) - {Ml + P 2 



01 Otl\/f2 i E>2\ 



Pchl, (10) 
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and hence Eq (Equation (|6|)). We then reconstructed the B hadron energy, E r B ec = 
E + E c h. Events were retained in which both reconstructed B energies satisfied E r B c < 
60 GeV. A final sample of 19,809 events was obtained with an estimated bb selection 
efficiency of 21.7% and a background contribution of only 0.23%, which was almost 
entirely from Z° — > cc events. 

The energy resolution of the final B sample is shown in Fig. 0, where we plot the 
normalised residual on Eb' (E B ue — E r B ec ) / 'E l B ue . The distribution was fitted to a dou- 
ble Gaussian function for which the mean positions, widths and normalisations were 
allowed to vary. 79.5% of the population lies in the 'core' Gaussian, of width 21.3%; the 
remaining population is characterised by a 'tail' Gaussian of width 31.3%. 

The angular resolution was similarly investigated and is shown in Fig. |3|, where we 
plot the residual on <p: (fr true — tp rec . The vast majority of angles are reconstructed to 
better than 5 mrad. The tail in the resolution function corresponds to B decays with 
shorter decay lengths, where the vertex position error had a larger relative effect on the 
determination of 0. 

The double-inclusive distribution of raw 5-hadron energies is shown in Figure 
Since we do not distinguish between B and B hadrons, the reconstructed S-decay ver- 
tices were labelled arbitrarily '1' and '2' on an event-by-event basis. 

We divided E B ec by the beam energy, E beam = y/s/2, to obtain the reconstructed 
scaled 5-hadron energy, x r B c = E B ec / Eb eam . 



4 Angle-Dependent B-B Energy Moments 

In each event we quantified the correlations between the two B hadrons in terms of the 
angle dependent scaled-energy moments proposed in Jl . We formed the single- inclusive 
I?-energy distribution and evaluated the moments (Eq. [I]) from the raw measured dis- 
tribution: 

d? c = ^M^TT 1 (ii) 

where the index k (1 < k < 20) runs over the bins of x r B c and iV is the number of events 
in the sample. Similarly, from the double-inclusive 5-energy distribution (Fig. (D we 
evaluated in each cos0 bin the double moments (Eq. 0): 

WW = ^^(x^y-^x^y- 1 (12) 

where the sum extends over the set of events in each cos0 bin. The normalised moments 
Gij (Eq. H) were evaluated: 

GTf(0) 

The first six moments, i = 1,2,3 and j - 
were defined by taking the average value 



Y)rec J~\rec 



(13) 



= 1,. . .,i are shown in Fig. ||. The bin centers 
of cosd> within each bin. 
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Also shown in Fig. |5| is a comparison with the simulated normalised moments; the 
simulation reproduces the data. Given that we showed previously |E|, that the Peterson 
function implemented in our simulation does not provide a good description of the 
6-quark fragmentation function, this agreement may naively appear to be surprising. 
However, if, as proposed in [[3], the non-perturbative contributions to the normalised 
quantity cancel, the agreement should be excellent, as observed. 

We used our simulated event sample to correct for the effects of the detector accep- 
tance, the efficiency of the technique for reconstructing 5-hadron decays, the energy 
resolution, and bin migrations caused by the finite angular resolution. We defined a 
binwise correction factor: 

Fij (0) = wmuc (14) 

where Gfj n ((p) \mc is the normalised moment calculated using generated true e + e~ — > BB 
+ X events, and G[ J ec (0)|MC is the corresponding moment calculated after simulation of 
the detector and application of the same analysis as applied to the data. For cos</> ~ — 1 
the value of is close to unity. As cos0 increases towards 1, rises monotonically to 
approximately 8 (-Fii) or 2 (-F33). The increase with cos0 reflects our decreasing efficiency 
to select BB events as the angle between the two S-decay vertices becomes smaller and 
the B energies decrease. For a given cos</> bin, F^ decreases as i and j increase due to 
the effective weighting of this efficiency by x^x^ 1 . 
We derived the true normalised moments: 

<W) = (15) 

which are shown in Fig [| The uncertainties associated with this correction procedure 
are discussed in the next section. 



5 Error Estimation 

We used our simulation to estimate the statistical error on the moments. The simulated 
event sample was divided into subsamples, each comprising the same number of events 
as in the data sample. The entire analysis was performed on each of these subsamples. 
Within each cos0 bin the r.m.s. deviation of the ensemble of results was evaluated and 
taken as the statistical error in that bin. For a given pair of i and j values the errors 
on the particular moment are not correlated between different bins of cos</>, but within 
a given cos0 bin the errors on the different moments are correlated, since the same set 
of events was used. 

We considered sources of systematic uncertainty that potentially affect our measure- 
ments. These may be divided into uncertainties in modeling the detector and uncer- 
tainties on experimental measurements serving as input parameters to the underlying 
physics modeling. For these studies our simulation was used. 
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Since our energy reconstruction technique is strongly dependent on charged-track 
properties, four sources of uncertainty were investigated: our simulated tracking effi- 
ciency, transverse momentum resolution, and the resolutions on the track impact pa- 
rameter and polar angles. In each case the simulation was corrected so as to reproduce 
the data. The full analysis was repeated on the data, and half the difference between 
the results obtained using the corrected and uncorrected simulations was taken as a 
symmetric systematic error. 

A large number of measured quantities relating to the production and decay of charm 
and bottom hadrons are used as input to our simulation. In bb events we considered the 
uncertainties on: the branching fraction for Z° — ► bb; the rates of production of B ± , B° 
and B° s mesons, and B baryons; the lifetimes of B mesons and baryons; and the average 
B-hadron decay charged multiplicity. In cc events we considered the uncertainties on: 
the branching fraction for Z Q — ► cc; the charmed hadron lifetimes, the charged multi- 
plicity of charmed hadron decays, the production of K° from charmed hadron decays, 
and the fraction of charmed hadron decays containing no 7i~°s. We also considered the 
rates of production of secondary bb and cc from gluon splitting. The uncertainty on the 
world-average value of each quantity was used to rederive the corrected Gy(0). In 
each bin of Gij(4>) the deviation between the rederived and standard values was taken 
as an estimate of the corresponding systematic error. Other relevant systematic effects 
such as variation of the event selection cuts and the assumed 5-hadron mass were found 
to be very small. 

The error associated with the choice of 6-quark fragmentation function used in the 
simulation was estimated by rederiving the corrected data using in turn the UCLA, 
Kartvelishvili and Bowler fragmentation functions 0. In each cos0 bin the r.m.s. de- 
viation of the results w.r.t. the standard value, using the Peterson function, was taken 
as an estimate of the systematic error. This error was typically much smaller than that 
arising from the other error sources. 

For each systematic error investigated, the shape of the angular dependence of the 
moments was not significantly changed, only the overall normalization was affected. 
For each moment, in each cos0 bin the dominant errors were typically those related 
to the uncertainties on the charged-track properties. The errors due to charm and 
bottom physics modeling were typically an order of magnitude smaller. In each cos</> 
bin all sources of systematic uncertainty were added in quadrature to obtain the total 
systematic error. The fully-corrected Gij, with statistical and systematic errors added 
in quadrature, are shown in Fig. || 

6 Comparison with Perturbative QCD Calculations 

The fully-corrected data were compared (Fig. ^ with a recent calculation M of the 
normalised moments complete at NLO in pQCD. Both the LO and full NLO calculations 
are shown in Fig. the calculations assume an a s (Mf) value of 0.120 and a pole b- 
quark mass value of 5.0 GeV/c 2 . It can be seen that the difference between the two is 
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relatively small, and that the LO calculation lies systematically slightly below the NLO 
calculation. The small size of the NLO relative to the LO contributions is an indication 
that the normalised moments are perturbatively robust observables. For each moment 
shown, the LO calculation undershoots the data. The NLO calculation reproduces the 
data across the full range of cos0. 

This comparison does not rely on any convolution of the pQCD calculations with 
models of the non-perturbative hadronisation process. Hence the excellent agreement 
between the pQCD calculations and the data verifies the ansatz of factorization be- 
tween the perturbative and non-perturbative phases that forms the basis for the pQCD 
calculation of heavy-hadron properties. 
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Figure 1: Pj-corrected mass (see text); data (points) compared with the simulation 
(histograms) in which the primary flavor content is indicated. 
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Figure 2: Resolution on the reconstructed 5-hadron energy. 
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Figure 3: Resolution on the reconstructed angle between the two B hadrons. 
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Figure 4: Reconstructed double-inclusive B-hadron energy distribution. 
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Figure 6: Corrected moments of the two Z?-hadron energies compared with LO and 
NLO pQCD calculations. The error bars are the sum in quadrature of the statistical 
and systematic errors (see text). 
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